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^~~^ We present updated predictions of effective elastic nucleon-nucleon cross sections 

^N intended for use in nucleus- nucleus reactions. A novel characteristic of the present 

^H approach combines all microscopic medium effects included in the Dirac-Brueckner- 
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Hartree-Fock G-matrix with a Pauli blocking mechanism which is more appropriate 
for applications in ion-ion reaction models as compared to a previous approach. The 
effective in-medium cross section is found to be quite sensitive to the description of 
Pauli blocking in the final configurations. 

I. INTRODUCTION 



> 

Q\ The investigation of the effective nucleon-nucleon (NN) interaction in dense hadronic 

O 

\0 matter is a topic of fundamental importance for nuclear reactions at intermediate energies 

^ (20 MeV/nucleon ^ Eiab ~ 300 MeV/nucleon) and for nuclear structure in general. The 

,__( relevant literature is very vast. Reference [1] is just a representative example of the tra- 

_^ ditional microscopic approach where two-nucleon correlations in nuclear systems are intro- 

duced through the G-matrix. Moreover, the effective NN interaction is the main ingredient 
of microscopic predictions of the nuclear equation of state (EoS) and thus impacts the prop- 
erties of compact stars. Dense hadronic matter can also be created in the laboratory in 
energetic heavy-ion (HI) collisions. Simulations of HI collisions are typically based on trans- 
port equations and describe the evolution of a non-equilibrium system of strongly interacting 
hadrons undergoing two-body collisions in the presence of a mean field. The Boltzmann- 
Uehling-Uhlenbeck equation ^ |3] and quantum molecular dynamics |1], along with their 
relativistic counterparts [5H7] , have been typically employed to describe intermediate-energy 
HI reactions. In-medium two-body cross sections are therefore an important component of 



such simulations. 

In direct reactions at intermediate energies the NN cross sections are often used as input 
to obtain quantum refractive and diffractive effects, replacing the role of optical potentials 
commonly used in low energy reactions [8j. Examples such as knockout (stripping and 
diffraction dissociation) reactions, elastic scattering, charge-exchange, and excitation of gi- 
ant resonances, are often carried out using reaction mechanisms based on the construction 
of scattering matrices built from the underlying NN scattering. Reaction calculations at 
intermediate to high energy are often conducted within the framework of the Glauber ap- 
proximation [H] and have been a frequent tool for testing nuclear models and constraining 
nuclear sizes. In fact, the description of complex nuclear reactions at intermediate energies 
based on individual NN collisions has a long tradition. In the framework of the Glauber 
model, the reaction cross section is written in terms of the "thickness function", which is 
the product of the averaged NN cross section and the overlap integral of the target and 
projectile local densities. 

In-medium NN cross sections have been calculated with a variety of methods. In semi- 
phenomenological approaches, one makes the assumption that the transition matrix in the 
medium is approximately the same as the one in vacuum and that medium effects come 
in only through the use of effective masses in the phase space factor [T0] - [T2] . Then, the 
in-medium cross section is scaled (relative to its value in vacuum) as the square of the ratio 
of the (reduced) masses. Phenomenological formulas, such as the one in Ref . |T3] , have been 
developed for practical purposes and combine the energy dependence of empirical free-space 
NN cross sections with the density dependence of some microscopic models. 

Microscopic predictions based on a medium-modified collision matrix were reported, for 
instance, in Ref. [14j, where Dirac-Brueckner-Hartree-Fock (DBHF) medium effects were ap- 
plied to obtain a medium-modified if -matrix. More recent microscopic calculations applied 
DBHF medium effects to produce a complex G-matrix including consideration of isospin 
dependence in asymmetric nuclear matter [T5] . 

It is the purpose of this paper to present our updated predictions of microscopic in- 
medium elastic NN cross sections with an improved description of Pauli blocking. The main 
objective is to produce two-body cross sections which include, microscopically, all important 
medium effects and are suitable for realistic applications in nucleus-nucleus scattering at 
intermediate energies including direct and central collisions. As explained in Section II, we 



start from a one-boson-exchange NN potential, which describes well the elastic part of the 
NN interaction up to high energy. Thus, as long as we are not interested in pion production, 
which is negligible up to, at least, several hundreds of MeV, it is reasonable to use NN 
elastic cross sections as input to the reaction model. Of course, the elastic part of the NN 
interaction can and does generate inelastic nucleus-nucleus scattering. 

In Section II, we describe the details of the calculation and highlight the differences with 
our previous approach. We then present a selection of results (Section III) followed by our 
conclusions and outlook (Section IV). 

II. DESCRIPTION OF THE CALCULATION 
A. The Dirac-Brueckner-Hartree-Fock G-matrix 

The starting point of our calculation is a realistic NN interaction which is applied in the 
nuclear medium without any additional free parameters. We use relativistic meson theory, 
which we find to be an appropriate framework to deal with the high momenta encountered in 
dense matter. In particular, the one-boson-exchange (OBE) model has proven very successful 
in describing NN data in free space and has a good theoretical foundation. 

The OBE potential is defined as a sum of one-particle-exchange amplitudes of certain 
bosons with given mass and coupling. In general, six non-strange bosons with masses below 
1 GeV/c^ are used. Thus, 

v= E KP^^ (1) 

a=7r, 77, p, a;, (5, cr 

with TT and t] pseudoscalar, a and 5 scalar, and p and uj vector particles. For more details, see 
Ref . Iin] . Among the many available OBE potentials, some being part of the "high-precision 
generation" [TTl HB], we seek a momentum-space potential developed within a relativistic 
scattering equation, such as the one obtained through the Thompson [19] three-dimensional 
reduction of the Bethe-Salpeter equation [20] . 

First, a self-consistent calculation of (symmetric or asymmetric) nuclear matter is per- 
formed within the DBHF approach [21]. This step yields, along with the EoS, the self- 
consistent nuclear matter potential, which is conveniently parametrized in terms of nucleon 
effective masses (see Ref. [H] for details). Then, the Thompson equation is solved for two nu- 
cleons scattering at some positive energy in the presence of a mean field due to the medium. 



The presence of the medium is accounted for through the (previously calculated) effective 
masses (applied in the two-nucleon propagator and also in the Dirac spinors representing the 
nucleons, consistent with the DBHF philosophy) and the presence of the (angle-averaged) 
Pauli operator to account for Pauli blocking of the intermediate states. 

In the usual free-space scattering scenario, the two-body cross section is typically rep- 
resented as a function of the incident laboratory energy, uniquely related to the nucleon 
momentum in the two-body center-of-mass frame, q, through relativistic invariants which 
yield the well-known relation Eiab = 2q^ /m. In nuclear matter, though, the Pauli operator 
depends also on the total momentum of the two nucleons in the nuclear matter rest frame. 
For simplicity, in the past we have used in- vacuum kinematics to define the total momentum 
of the two-nucleon system (that is, we assumed that the target nucleon is at rest, on the 
average). Schematically, the effect of Pauli principle on intermediate states arises in the 
G-matrix through the in- medium scattering equation [22] : 

. ip. ., . , ivi X [ ^V (q|v|qOQ(q^p)(q1G(p)|qo) ,^. 

(k|G(p)|qo) = (q|V|qo) - / ^^ E{pM)-E,-re ' ^'^ 

with qo, q, and q' the initial, final, and intermediate relative momenta of the NN pair in 
their center of mass, and p their total momentum. E is the energy of the two-nucleon system 
in the center-of-mass, and Eq is the same quantity on-shell. 

To account for Pauli blocking of the final state, we define the total elastic cross section 
as 

a-ivjv(g) = y (^) Q{q,p,9,p) dn , (3) 

where (da/dQ)^^^^ is the elastic differential cross section obtained from the G-matrix 
amplitudes as described above. 6 is the scattering angle and kp the Fermi momentum. The 
presence of the Pauli operator in Eq. ^ signifies that the integration domain is restricted 

by [15] 

kl-p^- q^ p^ + q^- kl 

< cos < . (4) 

2pq 2pq 

Setting Q = 1 in Eq. ([s]) amounts to ignoring Pauli blocking of the final state. (The virtual 
intermediate states are always subjected to Pauli blocking during the G-matrix calcula- 
tion which produces the amplitudes contained in (da/dQ)^^^^.) Additional simplifications 
result from the assumption that the differential cross section is isotropic. 




FIG. 1: Geometrical representation of Pauli blocking. 



B. The average in- medium cross section and Pauli blocking effects 



The cross section defined in Eq. ([s]) refers to an idealized scenario where a projectile 
nucleon, with some momentum above the bottom of the Fermi sea, strikes a target nucleon 
while both embedded in an infinite medium. For application to a realistic nucleus-nucleus 
scattering scenario, it is best to consider the situation depicted in Fig. [l] where the two 
Fermi spheres represent the local densities of the target and projectile ions, k is the incident 
momentum (the momentum of the colliding nuclei relative to each other), whereas ki and 
k2 are the momenta of any two nucleons. That is, ki and k2 + k are the momenta of the 
two nucleons with respect to the same point. Then, the relative momentum 2q and the total 
momentum 2p are given by 2q = k2 + k — ki, and 2p = ki + k2 + k, respectively. The larger 
circle in the figure is centered at p while |q| is the radius of the scattering sphere. The vector 
2q can rotate around the scattering sphere while maintaining constant magnitude due to 
energy-momentum conservation. 

Notice that, with the definitions given above, relative momenta which are off the symme- 
try axis of the two Fermi spheres (the k direction) are allowed, which is not the case with 
the assumptions made in Eq. (pi). That is, the two interacting nucleons can have momenta 
in arbitrary directions. In turn, this impacts the solid angle allowed by Pauli blocking, as 
shown below. (For completeness, we provide a detailed derivation of the allowed solid angle 
in the next subsection.) 
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In preparation for the nucleus-nucleus calculation, it is shown in Refs. [231 121] that the 
average NN cross section (assuming isotropy of the NN differential cross section) can be 
written as 

1 /■ 2(7 f 

^NN{k) = / rfki(ik2— o-AfAr(g) / dfl , (5) 

VfiVf2 J 1^ J Pauli 

where ki and ^2 are smaller kpi and kF2, respectively, and the angular integrations extend 
over all possible directions of ki and k2 allowed by Pauli blocking. Often, the empirical 
free-space NN cross section is used in the integral. In our case, <7nn{(i) = <J^^^^{q) is the 
(microscopic) NN cross section which contains additional medium effects as described in the 
previous subsection. Vpi and Vf2 are the volumes of the two (in general different) Fermi 
spheres. Because of azimuthal asymmetry, Eq. ^ can be reduced to a fivefold integration. 
Notice that the "symmetric" choice q = p = k/2 [25] amounts to making the approximations 
we adopted when writing Eq. ^. 

Finally, for an actual nucleus-nucleus scattering with given E/A, the average cross section 
given above becomes a function of the laboratory energy, E{k), and the local densities of the 
colliding nuclei, pi = 2kp^/{37c^), and are ready to be used in typical high-energy calculations. 
This is usually done by defining the average nucleon-nucleon cross section at the distance of 
closest approach b between the projectile and the target as 

/ fjphW Jd^n pi{ri)p2{ri + h) (Xnn{E,Pi,P2) .„. 

{crNN[E,b)) = — r — — — — — , (bj 

/<i3ripi(ri)p2(ri + b) 

where pi is the local density (at point r) inside nucleus i and aNN{E, pi, P2) is the in- medium 
NN cross section. 

The calculation of reaction cross sections in high-energy collisions is best described in the 
eikonal fomalism. The "survival amplitudes" (or S-matrices) in the eikonal approximation 
are given by [Qj [23] 

{aNN{E,b)) 



Si{E,b) = exp 



dq q pi (q) p2 (q) Jo (qb) 



(7) 



Air Jo 

where pi^2 (q) is the Fourier transform of the nuclear densities of the projectile and target, 
and the reaction cross sections are 

aR = 27r fdbb[l-\S{b)f]. (8) 



Applications to stable and unstable nuclei using Eqs. (l6jj8j) and our new prescription of 
Pauli-blo eking effects will be the subject of a future work. 




FIG. 2: (Color online) Two-dimensional projection of the geometry of Pauli blocking. 



C. Derivation of the Pauli-allowed solid angle 



As mentioned in the previous section, the relative momentum 2q and the total momentum 

2p are given as 2q = k2 + k — ki, and 2p = ki + k2 + k. We also define a vector 2b as 

2b = k2 + ki — k. Assuming that the collision is elastic, conservation of energy and 

momentum requires 

2p = k'l + k'a + k 

2q' = k2 - k'l + k (9) 

2b = k'l + k2 - k . 

The quantities kj^ and k2 are the momenta of two nucleons after the collision, whereas q 
is the relative momentum after collision, with |q | = |q|. Because of the Pauli exclusion 
principle, the following restrictions apply: 



or, 



|ki| = |p-q I > /cfi 
Ik'al = \h + ci\> kF2 , 



p^ + q^ — 2pq cos ai > kp^ 



(10) 



(11) 



&^ + g^ + 2bq cos ^2 > kp2 ■ 
In the equations above, ai is the angle between p and q , and ^2 the angle between b and 
q'. As illustrated in Fig. |2| we have 

q^-kl^ . b^ + q'^ 



cos 9/ 



p- 



2pq 



cos 9 



B 



Ivp2 



2bq 



(12) 




FIG. 3: (Color online) Pauli blocking of two nucleons in three dimensions. 



with 9a and 6b are the excluded polar angles. The excluded solid angles for each nucleon 
are then given by 

fia = 27r(l-coseA), fib = 27r(l -cos^b) , (13) 

and therefore the total allowed solid angle can be obtained from: 



n 



pauli 



A-K - 2{na + (^fc - fi) 



(14) 



where represents the intersection of the two conical sections Via and Vti,. The full calcula- 
tion has already been done in Ref. (2^; however, in this paper we will use a slightly different 
approach to calculate Vl. Figure [s] shows how Via and Vt^ are projected on the surface of a 
unit sphere. If Vti is the intersection of Via and i^f,, it is obvious that 



iie>eA + eB] 



a = i Qb iieB< Oa, e<\eB- ^ai ; 



na iieA<eB,9<\9B 



(15) 



The case \9b ~ 9a\ <9<9a + 9b''^^ more complex than the other three cases and a more 
detailed study is needed. As shown in Fig. Ill P and B are the centers of the two circular 




FIG. 4: (Color online) A different view of Pauli blocking of two nucleons in three dimension 



projections Qa and Qfy. The two circular contours intersect at R and L. a/2, /3/2 and 7 are 
the internal angles of the spherical triangle PBR. The circular sectors of Qa and Qb have 
areas equal to aQa/i^n) and f3Qb/{2iT), respectively. The intersection area of fia and Qb is 

given by 

a B 

Qi = —Qa + ^^b - 2ApnB . (16) 

Here, Ap^B is the area of the spherical triangle PBR. To obtain an expression for a/2, first 
we define the center of the unit sphere, O, as the orgin of the system, and Xp along the 
2;-axis. Point B is at location (l,6',a;/2), while point L has coordinates {1,6a,0). We can 
then write: 

OB ■ OL = cos 9b = cos 9a cos 9 + sin 9a sin 9 cos(a/2) , (17) 

from which a/2 can be readily obtained as 

' cos 9 b — cos 9 cos 9 a \ 



a/2 = arccos , 

y sm 9 sm 9 a 

In a similar fashion we find /3/2 to be given by 

cos 9 a — cos 9 cos 9 b 



/3/2 



arccos 



sin 9 sin 9b 
Applying the law of cosines of spherical trigonometry. 



(18) 



(19) 



cos 7 = — cos(q;/2) cos(/3/2) + sin(Q;/2) sin(/3/2) cos^ , 



(20) 
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FIG. 5: (Color online) In-medium pp cross section calculated as in Eq. (pi) for a variety of symmetric 
{kp\ = kp'i) and asymmetric (/ci?i 7^ kp'i) situations. 



we obtain 

7 = arccos(— cos(a/2) cos(/3/2) + sin(«/2) sin(/3/2) cos6') . 

From Girard's theorem of spherical trigonometry, we have 



(21) 



(22) 



Inserting Eq. (21 ) and Eq. (22 ) into Eq. ( 16 ), the solid angle Vti is found to have the following 
value 



fij = 2{7r-cos6'ACos ^{5ab)-co^0bco?, ^((5ba)-cos ^[co^9^{l - 5\b){1 - 5%j^-5ab^bA} 

(23) 

where 

cos 6i — cos 6 cos 6j 



Sij 



sin 9 sin 6j 



(24) 



Noticing that, while + Oa + Ob > t^, ^a and Vtb have two intersections on the hemisphere, 
we have 

= n,{e, Oa, Ob) + fi,(7r - 6, Oa, Ob) • (25) 
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FIG. 6: (Color online) As in the previous figure, for np scattering. 
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FIG. 7: (Color online) pp (left) and np (right) in- medium cross sections calculated from Eq. (pi) 
with kpi = kp2 = 1.0 fm~ . Solid red: predictions as in Figs.[5JJ6J dashed blue: the input NN cross 
section is evaluated in free space. See text for details. 



III. RESULTS 



A. Effective NN cross sections 



We begin by showing in Fig. |5] the average in-medium pp cross section calculated as in 
Eq. ([5]). On the left, we display a variety of cases with equal Fermi momenta, whereas 
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FIG. 8: (Color online) pp (left) and np (right) in-medium cross sections with kpi = kp2 = 1-0 fm 
Solid red: predictions as in Figs. [5]J6J dashed blue: predictions obtained with Eq. ([3]). 

asymmetric cases are shown on the right. Figure [6] contains the same information for the 
np cross section. After "overcoming" complete PauU blocking, the cross section generally 
rises with increasing incident momentum. In the np case, we observe, at least at the lower 
densities, a tendency to reach a broad maximum. In all cases, the cross sections become 
nearly flat at the larger momenta and begin to approach the free space predictions. 

Figures [5]J6] are more insightful when compared with Fig. [7j There, the pp and np cross 
sections shown by the dashed blue line are also calculated with Eq. ([s]), but the input 
NN cross sections in the integrand are evaluated in free space. Thus, comparing the two 
curves on the left-hand-side (or on the right, for np) shows the impact of the additional 
medium effects (besides those coming from the Jpaund^ factor in Eq. ([S])) originating from 
the G-matrix calculation and included in (Ttvat in the case of the solid curves. The impact 
is noticeable, with the microscopic medium effects further suppressing the cross section and 
shifting the position of the peak. We have chosen a particular case (kpi = kF2 = 1-0 fm~ ) 
for the purpose of demonstration, but the trend is similar for other densities. 

Figure [S] is crucial for the point that we wish to make. There, for pp and np (on the left 
and right, respectively), we compare the cross sections calculated from Eq. ([s]) to the corre- 
sponding ones evaluated with Eq. (pi) instead. The predictions from Eq. ([s]) have a sharper 
rise from zero and a more pronounced peak structure. As is reasonable, differences are large 
at low momenta, where the scattering is most sensitive to the description of Pauli blocking, 
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particularly near the onset of the cross section. Again, we have taken a representative case, 
but this pattern is common to all densities. It will be interesting to explore the impact of 
such differences on reaction cross sections, our next objective. 

IV. CONCLUSIONS 

Pauli blocking is perhaps the most important mechanism impacting the collision of two 
fermions in the medium. It is known to have a substantial effect on the scattering probability, 
that is, the in-medium cross section. 

In this paper, we predict in-medium effective NN cross sections suitable for applications 
to nucleus-nucleus scattering. The microscopic NN elastic cross sections, modified by all 
medium effects implied by the Dirac-Brueckner-Hartree-Fock theory of nuclear matter, are 
properly averaged so as to account for all possible directions of the relative momentum of two 
nucleons in the two colliding Fermi spheres. The more realistic description of the collision 
geometry amounts to an improved description of Pauli blocking as compared to a previous 
approach [15]. We find the effective NN cross sections to be very sensitive to the description 
of the Pauli blocking geometry. 

Our future plans include the application of these cross sections in Glauber reaction calcu- 
lations with stable and unstable nuclei. In closing, we also note that in-medium cross sections 
are related to the mean-free path of a nucleon, a fundamental quantity in the description of 
a nucleon propagation in nuclear matter. 
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